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Abstract. We study the dynamics of single-photon absorption by a single emitter 
coupled to a one-dimensional waveguide that simultaneously provides channels for 
^ ' spontaneous emission decay and a channel for the input photon. We have developed 

^ , a time-dependent theory that allows us to specify any input single-photon wavepacket 

^^' guided by the waveguide as initial condition, and calculate the excitation probability 

of the emitter, as well as the time evolution of the transmitted and reflected field. 
For single-photon wavepackets with a gaussian spectrum and temporal shape, we 
obtain analytical solutions for the dynamics of absorption, with maximum atomic 
excitation ~ 40%. We furthermore propose a terminated waveguide to aid the single- 
f^ ' photon absorption. We find that for an emitter placed at an optimal distance from 

Cn ■ the termination, the maximum atomic excitation due to an incident single-photon 

wavepacket can exceed 70%. This high value is a direct consequence of the high 
spontaneous emission /3-factor for emission into the waveguide. Finally, we have also 
explored whether waveguide dispersion could aid single-photon absorption by pulse 
shaping. For a gaussian input wavepacket, we find that the absorption efficiency can 
be improved by a further 4% by engineering the dispersion. Efficient single-photon 
absorption by a single emitter has potential applications in quantum communication 
and quantum computation. 
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1. Introduction 

Ultimate control over single light quanta, the emission of single photons, the absorption 
of single photons and the routing of photons between qubits, is of core interest for 
quantum information technology and extensively studied in [H El lU HI [5l E] . Ideally, 
one could use single photons that rarely interact with each other, and the environment 
as natural messengers of quantum information between nodes, where actual operations 
take place. Such nodes could then consist of localized atoms or quantum dots that can 
interact strongly with each other or external stimuli via, e.g., side band coupling, RF 
fields, or electrical gate signals, in an efficient and controllable way. Thus, it is extremely 
desirable to map a flying photonic qubit state onto an atomic qubit with unit probability 
[TJ [8]. However, such interfacing between an atom qubit and a flying photonic qubit is 
challenging, since it requires simultaneously an open photonic system for easy interfacing 
with freely propagating photons, yet also a high light-matter interaction strength, which 
is usually associated with the use of a high-Q, closed photonic system surrounding the 
atomic qubit. In absence of a high-Q cavity, one can either use a large ensemble of 
atoms to compensate for weak optical transition strengths of single atoms [9], or use 
highly focused optical beams to excite single atoms or molecules as efficiently as possible 
[TOj [TT| 112] . Among different mapping techniques, direct absorption of single photons 
by light emitters is an attractive option that may be realized using recent advances in 
the engineering of complex photonic environments [TOl [131 [3 [13 [HI [13 [HI [13 120] • 
Such efficient single-photon absorption is not only of relevance in quantum optics, but is 
also of central relevance for detectors, photovoltaics, or optical sensing and microscopy 
based on absorption or fluorescence [2T1 122] . 

Recent advances in nanophotonics enable the funneling of almost all the single 
photons emitted by a single emitter into a single mode [231 (211, or directing single 
photons into narrow beams [251 ESI ISZj and even dynamical steering of singe-photon 
emission |28]. As sketched in Fig. [H it has been shown that these one-dimensional (ID) 
or quasi- ID waveguides can be used to efficiently control the spontaneous emission 
[23] [2H l29l l30l [3T1 [32] . Inversely, reciprocity in classical electrodynamics predicts 
efficient coupling of single photons with an emitter, provided the incoming single-photon 
wavepacket is fed through the channel which is associated with a high spontaneous 
emission rate. Recent theoretical work [T71 [191 IMl IMl [35] has hence focused on the 
interaction of single emitters and single photons guided by ID waveguides. However, 
this work mainly concerns the reflection and the transmission probability of single- 
photon wavepackets of very narrow spectral bandwidth interacting with a single emitter. 
These models consist of stationary solutions for the interaction of a two-level system 
with a one dimensional waveguide. In such stationary cases, the photon wavepackets 
have a quasi- inflnite temporal extent, and hence the atom is essentially in the ground 
state at all times. In this case, the atom can simply be treated as a point scatterer 
[ini [111 [36] . To optimize the atomic excitation probability, a time dependent treatment 
is required to predict which photon wave packets have optimal temporal and spatial 
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Figure 1. Possible experimental realizations of a ID waveguide strongly coupled 
to a single emitter: (a) a propagating surface plasmon polariton (SPP) mode on a 
metallic nanowire, (b) a guided mode in a photonic crystal waveguide and (c) a guided 
plasmonic mode on an infinitely long plasmon particle chain. 



profiles for atomic excitation. Originally it was pointed out by Cirac et al. [SI EZ] that a 
time-reversal symmetric photon wavepacket can be used to efficiently transfer quantum 
states among distant nodes consisting of a A-type atomic medium. The interesting 
concept of time-reversal symmetry was also applied to two-level atoms, showing that 
it is indeed possible to perfectly invert an atomic qubit using photon wavepackets that 
are spatially and temporally the inverse of the photon wavepacket emitted by a qubit 
through spontaneous emission [201 EH] • However, realizing this prediction requires highly 
non-trivial pulse shaping, especially at the single-photon level [39]. Thus there is large 
need for a time-dependent theory that quantifies what atomic excitation efficiencies 
can be reached using practically achievable single photon wavepacket and photonic 
structures. 

In this paper we develop a time- dependent theory that quantifies what atomic 
efficiencies can be reached using practically achievable single photon wavepackets in ID 
waveguides. As sketched in Fig. (H we consider a single-photon wavepacket propagating 
along a ID waveguide, interacting with a single emitter. We explore the possibility of 
maximizing the single-photon absorption by the emitter via engineering the photonic 
environment of the emitter, and by shaping the pulse of the input single-photon 
wavepacket through waveguide dispersion. We study ID or quasi-lD waveguides, 
considering all the spontaneous emission channels, since it is the competition between 
the pumping channel and the spontaneous emission channel of the emitter that will 
ultimately determine the absorption efficiency. We focus on the role of the temporal 
coherence of the input single-photon wavepacket for exciting the emitter in the ID 
waveguide. The formalism can be generalized to treat a three-dimensional (3D) light- 
scattering problem for examining the role of spatial coherence, which will be addressed 
in a future paper. 

The paper is organized as follows. In Section 2, we formulate the model and derive 
general solutions for the dynamics. Section 3 presents results obtained for gaussian 
input single-photon wavepackets propagating along an infinitely long ID waveguide. We 
obtain an analytical solution for the dynamics of the atomic excitation. Furthermore, 
we present a physically transparent model that relates our dynamical model to the 
stationary reflection and transmission spectrum, and we present a simple model for the 
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time-dependent reflection and transmission. We find that the maximum absorption 
probabihty strongly depends on the temporal coherence of the input single-photon 
wavepacket. In Sec. 4, we consider an emitter located near the termination of a half- 
truncated ID waveguide. Near such a termination, the spontaneous emission /3-factor 
for uni-directional emission can be strongly enhanced. By reciprocity, uni-directional 
emission with a high beta factor also results in the highest possible atomic excitation 
probability. Our work thus provides a simple route to optimize photonic waveguide 
structures for photon-qubit interaction. Finally, we investigate the possibilities of 
improving the atomic excitation by pulse chirping, or pulse dispersion to modify the 
temporal coherence. Section 5 concludes the paper. 

2. Model 

In this section, we will outline our theory for calculating the dynamics of single-photon 
absorption. Our model follows a similar procedure employed by Dorner et al. [IQ] 
for spontaneous emission and resonance fluorescence of a single emitter coupled to a 
mirror. In section 12. ![ we construct a general model for single-photon absorption valid 
for arbitrary (non-absorbing) 3D photonic structures with arbitrarily shaped input 
wavepackets. In Section 12. 2[ we apply the general model to study the dynamics of 
single-photon absorption by a single atom coupled to a ID waveguide [HI [33]. The 
reflection and transmission coefficients derived in stationary theories can be extracted 
by evaluating our time-dependent solutions at times far later than the time interval 
where the photon wavepacket collides with the atom, and after the emitter has relaxed 
back into the ground state due to spontaneous emission. 

2.1. Hamiltonian, initial conditions and equations of motion 

In a finite volume the electromagnetic field can be decomposed into discrete modes. 
In the discrete-mode quantization scheme, the interaction between a single-photon 
wavepacket and the emitter is modeled by the Hamiltonian 

-^tot = Hq + Hj, (1) 

where ^o defines the free constituents Hq = huo\e){e\ + ^ /i(jJAa"''(A)a(A), and where 

A 

Hj defines the interaction part Hj = J2H9x^W\^){9\ ~^ 9x^''W\9){^\]- Here \g) and 

A 

|e) represent the atom ground and excited state respectively, uq denotes the atomic 
transition frequency, while a(A) is the annihilation operator for the mode A. Without 
loss of generality, we take the emitter position as our real-space origin (emitter located 
at r = 0). We expand the Schrodinger- picture state {ipit)) at time t in the basis of all 
states with one excitation, i.e., |e, 0), \g, 1a), as follows: 

\m)=Qit)e-'-'%0) + J2e-'-^'Cl{t)\g,l,). (2) 
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By separating the time dependence driven by Hq, the coefficients Co(^) ^^^ ^xi^) ^^^ 
essentially solved in the interaction picture. Since we are interested in the absorption of 
a single-photon wavepacket, we assume that in the initial state, the atom is unexcited 
(Cq(0) = 0), while the initial wavepacket is described by Cf (0) = ^(A). 
The equations of motion are 

^^dC^^ = J2Clit){e,0\Hj\g,We-^^-^-^'^\ (3) 

A 

'^^t = Co^(t)(^,l,|i//|e,0)e^(^-^«)*. (4) 

The transition matrix element {g, lA|-f^/|e, 0) can be simplified as {g, l\\Hi\e,0) = hg\, 
where g\ is the well-known interaction strength of a single emitter with mode A |43j . 
Integrating Eq. (jlj) yields 

Ciit) = CliO) - z /* C^,it')gle^^-^--'>^Ut\ (5) 

Jo 

and using this in Eq. (^, one can obtain 

A A "^0 

Equation (jG]) is of central importance in this work: solving this equation of motion 
provides the time dependent excitation amplitude of the single emitter, hence 
quantifying the efficiency of single-photon absorption. The right-hand-side (rhs) of 
Eq. (P) consists of two terms. The ffist term accounts for excitation of the atom by 
the single-photon wavepacket. The second term accounts for de-excitation of the two- 
level system by emission of a photon. Two facts are immediately obvious: ffistly, 
Eq. (|6]) reduces to the Weisskopf-Wigner theory of spontaneous emission when the 
atom is excited and all the modes of the field are empty initially |12]. Secondly, the 
maximum excitation probability of the atom depends on the shape and duration of 
the incident photon wavepacket. Indeed, the incident photon wavepacket is the only 
channel for driving the atom, while any excitation is continuously subject to exponential 
decay due to vacuum fluctuations. In the Weisskopf-Wigner approximation, the second 
term of the rhs in Eq. ([6]) is memoryless, and can be simplified to rCo(t)/2, where 
r = 27r ^ g\gx5{uj\ — Uq) is the total spontaneous emission decay rate. By decoupling 

A 

Co(t) from the time dependence of C^(t), namely, by combining the two ffist-order 
differential equations defined by Eq. ([3]) into a second-order differential equation defined 
by Eq. ([6]), one only needs to know Cf (0) in order to solve for Co(t). In Eq. ([6]), C^(0) 
represents the mode distribution in the incident single-photon wavepacket. The entire 
spatial structure of the photonic modes surrounding the emitter is implicitly encoded in 
the coupling strength gx^ which can be further expressed as g\ = - — ^'^ , where fj, is 
the dipole moment of the emitter, and where Exi{r = 0) is the normalized single-photon 
field strength for mode A, at the position r = of the emitter. 

In the above, we have assumed that the fields belong to a finite quantization volume. 
In a ID or quasi- ID waveguide system, a continuous- mode quantization scheme needs 
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to be adopted. We assume that the single-photon wavepacket supported by the ID 
waveguide can be decomposed into cyhndrical waves, i.e., E\{r) = Ex{x,y)e^^^^~'^^\ 
with the mode label A = {m,p, /9,g}, where (3 is the component of the wave vector 
along the z axis, q represents the magnitude of the wave vector perpendicular to the 
z axis, m is the angular momentum, and the index p is used to distinguish between 
two degenerate polarization modes for given m, /3, and q. For more details about 
the normalization of these continuous modes and the substitution of summations over 
discrete modes by integration over continuous modes see Appendix A. We also assume 
that the incoming single-photon wavepacket has a very narrow bandwidth, on the order 
of 10^ rad/s, since its temporal duration is comparable to the SE lifetime. Such a narrow 
bandwidth with respect to the atomic transition frequency ujq = 10^^ rad/s ensures the 
the vahdity of the linear dispersion relation, i.e., /S^^ — /3o = (w — ijJo)/vg, where Vg is 
the group velocity of the propagating mode. From Eq. ( lA.ll) and the Weisskopf-Wigner 
approximation, one obtains, 

^^ = --^M -[jdujY^ xlA^)E^,nxUr = 0)e-^-*] - rCo^(t)/2, (7) 



where xt ni^) — 7 ; ^^^ ^ i^ given by k = {m,p, g}. Equation (j7]) can be used to 
calculate the excitation probability of the single atom as it is illuminated by a single- 
photon wavepacket in a complex photonic environment, e.g., in a photonic cavity or 
near an optical antenna. 

2.2. Dynamics of absorption, reflection, and transmission in ID and quasi-lD 
waveguides 

In this section we consider the dynamic response of both the emitter and the photon field, 
taking as initial condition a single-photon wavepacket incident from the -z direction, 
traveling towards the emitter at 2; = 0. The transverse distribution of the waveguide 
mode is independent of z, and the longitudinal wave number of the modes k is assumed 
to obey the linear dispersion relation k — ko = {u — u^j/vg in the frequency range 
of interest. We study the dynamics of coherent single-photon absorption with special 
focus on the dependence of absorption on the properties of the input single-photon 
wavepacket. The dependent variable time t has a one-to-one correspondence with the 
position of the peak of the single- photon wavepacket. Initially, the peak of the pulse is 
at the position z = —Zq, and the interaction between the single-photon wavepacket with 
the emitter reaches its maximum at t = Z^/vg. At approximately this time, the emitter 
reaches its maximum excitation probability. Meanwhile, the emitter also loses atomic 
excitation probability due to the fact that it decays via spontaneous emission in both 
the forward and backward direction. The emitted light can interfere with the incident 
beam, which results in extinction and pulse reshaping in the transmitted channel. The 
relevant timescale for de-excitation by spontaneous emission is the lifetime l/F, where 
F is the spontaneous emission decay rate of the emitter. The single-photon wavepacket 
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will be redistributed into a (re-shaped) reflected and transmitted pulses and the emitter 
will relax to the ground state in several lifetimes after attaining the maximum excitation 
probability. 

In any ID waveguide, the continuous-mode variable can equivalently be taken 
as the wavevector k or the frequency Uk due to the linear dispersion approximation. 
The wavevector k is chosen here. In such infinitely long single mode waveguide, the 
spontaneous emission divides equally between the forward and backward propagating 
channel. If we assume that the incident photon wavepacket is a single packet offered 
from just one direction, this means that the spontaneous emission /3o-factor for the 
pumping channel is at most 50%. For any given incident photon wavepacket C^(0), it is 
straightforward to calculate the dynamic atomic excitation CQ{t) in the ID waveguide, 
which is given by, 

oo 



—I 



I dA;Cf (0)(7fce-*("'=-"°)* - rCo'(t)/2. (8) 



dt 

In our case, the coupling strength g^ is assumed to be frequency-independent due to the 



narrow bandwidth of the input single-photon pulse and given as gk = \ r^Vg [19]. 

Once the atomic excitation Cq (t) is known, one can calculate the dynamics of Cf,{t) 
to obtain the reflected and transmitted wavepacket. According to Eq. ([5]), the total 
probability amphtude of the forward propagating wavepacket is a coherent superposition 
of the initial incident photon wavepacket and the emitted photon, and is given by 

CI At) = cm - . f C^,{t')gle^^-'^-^^' dt' . (9) 



For the backward direction, there is no incident term (C|_(0) = 0). Hence, the 
corresponding probability amplitude is 

Ciit) = -I f C^,{t')gle^^--'-"^' dt' . (10) 



In most realistic cases, an emitter coupled to a ID waveguide will still have a residual 
coupling to free space modes that are not guided by the wire. For such a quasi- 
ID waveguide, the spontaneous emission /3-factor into the forward plus backward 
propagating waveguide modes will have a value less than 100%. Also, a break in 
symmetry in the geometry might imply that emission into forward and backward 
propagating waveguide modes is unbalanced. In these cases, where the pumping channel 
funneling the incident single-photon wavepacket has a spontaneous emission /3-factor of 
/So, Eq- dHl) becomes. 



dC-.Jt) _ ^ /■^,^^5.n^.//3o^ 
dt 



-I J rffcC,^(0)^^^,e-(— »)* - TC^,{t)/2. (11) 



It should be remarked that the first term of the rhs in Eq. f fTTj) is simply the pulse 
shape in absence of light-matter coupling, assuming that the coupling strength gk is a 
constant. In the next section we will study specific wavepackets and make use of the 
equation of motion derived here. 
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Figure 2. (a) Snapshots of the amphtudes of the transmitted and reflected 
wavepackets as distributed along the spatial axis at different instances of time, (b) 
the time dependence of the atomic excitation. At f = 0, the peak of the gaussian input 
single-photon wavepacket is localized at z = —Zq = —lOvg/T, with FWHM equal to 
r. The red dot indicates the fixed position of the emitter. 



3. Gaussian pulse excitation of one-dimensional waveguide 

First we consider gaussian input single-photon wavepackets defined as C^(0) 



at t 



46 



'i{k-ko)Zo- 



(fco-fc) 



where 

r-oo 



Zq is the position at which the peak of the pulse passes 
0. The relation J |C^(t = 0)\'^dk = 1 ensures that there is only a single-photon 
in the wavepacket. With the linear dispersion relation and the substitutions of to = ~ 
and fi = Afg, one can obtain an analytical solution for the dynamics of single-photon 
absorption, i.e., for the excitation amplitude of the two level system 



1 
'2 



—e " 
a 



b b 

'[erf(Vat + —^) - erf(^)]e" 
'a \ a 



(12) 



4 ' 



where s = —i\/-^[27rDI 



= -(r - nHo)/A, c = nHl/A, and erf(a:) is the 

error function, defined as erf(a;) = -y^ Jn g~^ dt. The result CQ(t) obtained from the 
time-dependent theory [1U| HT] contains much richer information than can be obtained 
from stationary theory, due to the fact that the initial conditions can not be included 
in the Fourier transform descriptions of light scattering by emitters [HH [191 SOI l4T] . 
However, we would like to point out that an alternative route to solve the problem that 
includes the initial condition is the Laplace transform method developed by one of us 
(See Wubs et al. |11]). We find that the two methods are equivalent and yield the same 
result for |CQ(t)| . The derivation of |CQ(t)p based on the Laplace transform method is 
briefly outlined in Appendix B. 

Figure |2t^a) shows the spatial variations of the forward and backward propagating 
wavepackets, which are obtained by Fourier-transforming C|_^(t) and Cf. _{t) with 
respect to k. Initially the incident single-photon wavepacket is far away {z < —Zq) 
from the emitter. At this instance there is negligible interaction. As the leading edge 
of the pulse reaches the emitter, the interaction is effectively switched on. From this 
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time onwards, the emitter gains excitation amplitude, and starts to emit an outgoing 
wavepacket into both the forward and backward directions. Before the emitter reaches 
its maximum excitation shown in Fig. EJ^b), one notices that the forward-propagating 
wavepacket experiences a sharp drop at the position of the emitter, while the backward- 
propagating wavepacket experiences a steady increase of its magnitude at z = 0. 
Figure E^b) shows the time dependence of the corresponding atomic excitation, which 
will be discussed in sub-section {3.2). The sharp drop in the forward propagating 
wavepacket is due to the transfer of energy from light to the emitter, the rate of which is 
larger than the SE decay rate of the emitter itself, given our assumption on the coupling 
strength gk, set by the mode profile and Vg. In contrast to the forward propagating 
packet, the backward propagating wavepacket obtained purely from emission has an 
increasing magnitude. After reaching the maximum excitation, the relaxing emitter 
leads to decreasing magnitudes for both forward and backward propagating wavepackets. 
Moreover, the forward-propagating wavepacket is seen to be strongly reshaped due to 
interference of the incident packet with the emitted light. The resulting minimum in 
the forward propagating wavepacket will be further discussed in the stationary limit in 
sub-section {3.1). One also notes that the trailing edges of the two wavepackets are 
longer than the leading edges. The widths of the leading edges are determined by the 
incoming single-photon wavepacket, while the trailing tails are essentially determined 
by the lifetime of the excited state of the emitter. 

3.1. Stationary limits of the single-photon absorption in ID waveguide 

We now examine our theory in the stationary limit, to connect our dynamic results 
to earlier work by Shen and Fan [33j- The stationary limit essentially corresponds to 
single-photon pulses with large temporal widths and low amplitudes, implying also that 
the emitter is always in the ground state. From our dynamic theory, we can extract 
transmission and reflection amplitudes in the limit of long excitation pulses, and at 
times far later than the time window in which interaction with the atom has taken place. 
Using a Fourier transform, we extract frequency dependent reflection and transmission 
amplitudes. We compare these frequency dependent coefficients, which we refer to as 
the stationary limit of our dynamic theory, with the frequency- dependent coefficients 
obtained in purely stationary theory by Shen and Fan p3]. In the stationary limit, 
in which the emitter has relaxed into its ground state, the forward- and backward- 
propagating wavepackets are essentially the transmitted and reflected light, similar to 
light scattering by an impurity in a ID waveguide. If we take a snapshot long after the 
interaction time (we use tg = lOO/F, but our results are independent of this choice), 
we find instantaneous reflection and transmission amplitudes in fc-space. By taking 
advantage of the linear dispersion relation, we can convert these A;-space instantaneous 
amplitudes into spectral information. The transmission and reflection spectra can be 
extracted as 

rM = |c^^,+(g/C4(o)p, (13) 
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Figure 3. (a) Transmitted spectral intensity and (b) reflected spectral intensity of 
a gaussian incident single-photon wavepacket mediated by a single emitter. Plotted 
results are for different ratios of pulse bandwidth fl to decay rate F. Note that the x-axis 
represents frequency in units of the pulse bandwidth in (a), (b) but in units of decay 
rate in (c). In (a), (b) we find that as the ratio of bandwidth to decay rate is decreased 
from large (5r, red), to small (0.2r, black), the reflected intensity appears to increase 
strongly. The transmitted intensity shows a pronounced minimum at zero detuning. 
This minimum is deepest, reaching T = 0, when the incident pulse width exactly 
matches the lifetime of the emitter (dark green, Jl = F). (c) shows the transmission 
and reflection spectra, obtained by normalizing the transmitted and reflected intensity 
as plotted in (a), (b) to the intensity spectrum of the incident wavepacket. The stars 
and triangles correspond to the results extracted from our time dependent formalism, 
via Eq. (12). The solid lines show the stationary calculations. 



^M = FL-(ts)/c^^,(o)r. (14) 

At first sight, the coexistence of the time and frequency dependence in the transmission 
and reflection spectra may appear odd. However, for times tg much later than the 
interaction time, the time tg is essentially equivalent to a spatial variable z = Vgtg- 
Hence, the transmission and reflection spectrum can be interpreted as frequency signals 
monitored at a certain position. 

Figure [3] shows the spectral intensities of the transmitted wavepacket and the 
reflected wavepacket, as well as the transmission and reflection spectra in the stationary 
limit, at ts, i.e., long after the interaction of the emitter with the incident wavepacket. 
Figure [3t^a) and (b) show that the spectral intensity of transmitted light has a dip at 
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the resonance frequency, while the spectral intensity of the reflected light resembles the 
original pulse shape. The dip corresponds to the minimum in the forward propagating 
wavepacket in Fig. |2]^a) due to the resonant interaction with the emitter, and the 
dip magnitude depends on the bandwidth of the incident wavepacket. We also note 
that narrower bandwidth will yield more reflection of light. Figure [^c) shows the 
transmission and reflection spectra in the stationary limit, by normalizing the reflected 
and transmitted pulse spectra to the spectrum of the incident wavepacket. We now 
compare the transmission and reflection spectra from our limiting procedure with 
stationary solutions obtained by solving for the eigenstates of the system as reported 
by Shen and Fan [33] . Such stationary behavior of the transmission and reflection can 
be well modeled by solving for the eigenstates of the system. For a given detuning 
6 = Uk — oJo, the corresponding reflection and transmission spectrum can be described 
as follows 1331. 



T{6) 



RiS) 



6 



iT/2 + 6 
iT/2 



(15) 



iT/2 + 6 ' 

The transmission and reflection spectra that our time- dependent theory predicts using 
a gaussian input wavepacket with narrow bandwidth agree well with Shen's stationary 
model deflned by Eq. ( !T5|) . as shown in Fig. [3] (c). We also like to point out that the 
transmission and reflection spectra for narrow bandwidth pulse in our stationary limit 
are independent of Q. This independence, which is rigorous only in the limit of zero 
bandwidth, is confirmed by the identical transmission and reflection spectra for the 
three pulses in Fig. |3] (a) and (b). Based on the stationary model, the full width half 
maximum (FWHM) of transmission and reflection equals the total decay rate of the 
emitter, i.e., FWHM=r. 

3.2. Time dependence of the transmittance and reflectance of the single-photon 
wavepacket in ID waveguide 

The stationary model accounts for the transmission and reflection at times long after 
the interaction for any given frequency distribution. Indeed, a good approximation of 
the transmitted and reflected wavepackets are obtained simply by multiplication of the 
incident spectrum, with transmission and reflection coefficients combined with the right 
phase factor. However, the stationary model does not account for the atomic excitation 
probability of the emitter during the interaction process. This dynamic information is 
an essential result of our model, that we now proceed to discuss. Also, the model allows 
to obtain the time dependent transmission and reflection probabilities at the time of 
interaction, instead of being limited to times much later than the interaction interval. 
We deflne the time dependent transmission (reflection) probabilities by including all 
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Figure 4. (a) Time dependence of the transmission and reflection probability of the 
incident single photon wavepacket and the atomic excitation, (b) Contour plot of the 
atomic excitation probability versus time and bandwidth of the incoming gaussian 
wavepacket. 



amplitude traveling to the right (left) as follows: 



T(t) 



rffc|Cf,(t)p 



(16) 



R(t) 



dk\Cl_it)\'. 



(17) 



It is important to realize here, that these coefficients are defined not via spatially 
separating the amplitude present to the left/right of the atom, but strictly by separating 
forward and backward direction by wave vector. As opposed to the refiection and 
transmission coefficients in the stationary limit, these instantaneous coefficients have 
no frequency content. Instead they allow us to monitor the distribution of excitation 
between the atom, and the forward and backward emission channel as a function of time. 
The frequency dependent transmission and reflection reported in Fig. |3] can be obtained 
by evaluating T(t) and R(t) for large times (e.g., t > lOO/F), for many narrow-band 
initial conditions centered at different frequencies. 

In Fig. m we consider the time-dependent probabilities for the atomic excitation, 
transmission and reflection. For an incoming single-photon wavepacket with FWHM 
equal to F, we flnd that the atomic excitation increases as a function of time as the 
peak of the pulse approaches . After t = 10/F, the atomic excitation decays with a time 
constant give by the lifetime of the emitter (1/F), after achieving a maximum value of 
atomic excitation close to 40%, see Fig. H] (a). The reflection builds up with time due 
to the continuous re-emission of light from the emitter. In contrast, the forward packet 
has a non-monotonic, though generally decreasing behavior. The kink in T{t) is due to 
interference of the incident packet and re-emitted light. In contrast to the prediction of 
the stationary model, the transmitted wavepacket does not vanish, although the incident 
wavepacket is tuned exactly to resonance. The imperfect reflection as compared to the 
perfect reflection in the stationary limit is due to the fact that optimum atomic excitation 
requires a flnite pulse length, or equivalently flnite bandwidth of the incident packet. 
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whereas perfect reflection only occurs at a single frequency. The impact of the finite 
bandwidth of the incoming wavepacket on the maximum value of the atomic excitation 
is also evident in Fig. |l](b), which shows the atomic excitation probability as function of 
time for different bandwidths of the incoming wavepacket. The graph clearly shows that 
there is a range of optimum bandwidths, that lead to atomic excitation probability with 
a maximum value of around 40%. This range of optimum bandwidths is comparable to 
r, corresponding to incident photon wavepackets that have a duration comparable to the 
spontaneous emission decay time. The 40% atomic excitation is surprisingly high given 
that reciprocity sets the fundamental limit for excitation from just one direction in the 
wave guide to be 50%. The value of 40% could be pushed closer to the limit of 50% by 
not using Gaussian pulses, but rather the inverting pulses proposed by Stobihska et al 
[H EZl EHl [SH]. However, exceeding the limit of 50% will invariably require illuminating 
the emitter from two sides with a proper phase relation between the two input pulses. 
In the following we show that instead of increasing the complexity of illumination by 
using two pulses, it is also possible to further increase the excitation probability above 
50% by engineering the photonic environment of the emitter to have broken symmetry. 

4. Gaussian input on a semi-infinite one-dimensional waveguide 

For an infinitely long ID waveguide, the atomic excitation probability is prevented from 
reaching unity by the fact that the emitter can decay equally into two directions, namely 
the forward and backward directions. Essentially, the coupling efficiency of the emitter 
with the one-sided pumping channel, or conversely the spontaneous emission /3-factor in 
just one waveguide direction, is at best 50%. In order to obtain higher atomic excitation 
probability, all the channels into which the emitter decays should be suppressed, except 
for the one optical input/output channel through which the incident photon wavepacket 
is sent in. In this section, we consider a situation where the single-photon wavepacket 
is launched into a semi-infinite ID waveguide. By placing the emitter at an optimized 
distance from the waveguide termination, we expect that the coupling of the emitter to 
the input waveguide mode is optimized. 

4.1. Enhanced spontaneous emission and the spontaneous emission (3-factor by a 
terminated metallic nanowire 

Specifically, we study a terminated metallic nanowire coupled to an emitter, as shown in 
Fig. [5](a). On the right-hand side, the metallic nanowire is terminated with a spherical 
endcap of radius R (equal to the radius of the nanowire). On the left-hand side, the 
nanowire is infinitely long. As shown in Fig. [5] (a), a is the distance of the emitter to the 
right rod end, and d is the distance to the wire edge. Figure [5t^b) shows the variation 
of the decay rate as a function of a for d=lQ nm, and an emission wavelength of 1000 
nm. We obtained numerical results using finite-element method (FEM) calculations, 
as explained in our previous work [32]. Such FEM modeling is rather fiexible and can 
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Figure 5. (a) Sketch of the terminated metalhc nanowire coupled to a singe emitter, 
(b) Position dependence of the total decay rate of the emitter coupled to the half 
truncated metallic nanowire. Here a denotes the distance of the emitter to the rod 
end on the right. The 20 nm radius silver nanowire (e^g=-50+0.6j) is embedded in 
the host material with index n—lAlA. The inset in (b) shows the total decay rate as a 
function of length of the modeling domain, which does not give any significant impact 
on r, more details see Ref. |32]. (c) Time dependence of the atomic excitation for a 
terminated metallic nanowire with spontaneous emission /3- factor of 91%. (d) Contour 
plot of the maximum of the atomic excitation versus bandwidth and /Ss for a gaussian 
input wavepacket in a terminated nanowire with /32 = 0. 



handle complex photonic structures |46] compatible with current lithographic fabrication 
technology. For an infinitely long 20 nm radius silver nanowire, there are only two 
guided plasmonic modes, i.e., a forward and a backward propagating guided mode with 
propagation constants ±k and corresponding wavelength Aeff = 27r/k =436 nm. The 
decay rate into the forward and backward propagating plasmon mode is 33.8, when 
normalized to the decay rate of the emitter in vacuum. Accordingly, the decay rate into 
one direction is 'jpi^Q ~ 17. The total decay rate coupled to an infinitely long nanowire, 
including the decay into radiative modes in the continuum, and direct quenching, is 
r=38.27. Hence the beta-factor into the pumping channel is approximately 44%. For 
a terminated metallic nanowire, the total decay rate shows a clear interference pattern 
as a function of the distance of the emitter to the termination. This pattern is exactly 
the ID analog of of Drexhage's observation that the lifetime of a molecule in front of 
a mirror oscillates as a function of mirror-emitter separation |45]. Compared to the 
infinitely long wire with r=38.27, we can enhance or inhibit F by a large amount, up 
to ~ 80%, by tuning the source-termination distance a. Importantly, by integrating the 
power flux in the transverse plane of the plasmon nanowire for a = 190 nm, we find that 
91% of the spontaneous emission is coupled into the guided mode that exits the plasmon 
wire on the left hand side. We hence also expect up to 91% excitation efficiency when 
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offering a single photon wavepacket through the waveguide. 

The mechanism responsible for the enhancement is that the waveguide termination 
essentially functions as a mirror, which reflects the right propagating plasmonic mode 
to the left. Similar to image dipole theory for Drexhage's experiment, which predicts 
a two-fold enhancement of the spontaneous emission decay rate for a dipole emitter in 
front of and perpendicular to a mirror, the waveguide termination also acts as a mirror, 
which yields an enhancement factor of 2 of the spontaneous emission into the plasmonic 
modes. Consequently, the emission rate into just the left propagating plasmonic mode is 
approximately enhanced by a factor of 4 since the two-fold enhancement is distributed 
over just half the channels. Compared to the infinitely long nanowire, such termination 
with a proper distance to the emitter gives rise to a spontaneous emission /3-factor up 
to 91% for the pump channel. 

4-2. Single photon wavepacket propagation along the terminated metallic nanowire 

Having quantified the spontaneous emission decay rate enhancement at the wire 
termination using FEM modeling, we now proceed to model the probability of absorption 
of a single photon wavepacket. To this end we first need to estimate the shape of 
the wavepacket including the reflected part, for which we need to know the complex 
reflection coefficient at the wire end. We can estimate the complex reflection coefficient 
from the calculated decay rate. Indeed, if the reflection coefficient at the rod end is 
re-'^, the decay rate of the emitter influenced by the rod end can be approximated 
as ^^ oc l^tlVl^oP = |^o(l + re^'^e^'2''»'')|Vl^oP, where Eq is the electric fleld for 
the plasmonic mode without reflections, Et is the total field including the left-hand 
propagating plasmonic mode and the one reflected from the right. From Fig. [5] (b), we 
extract the reflection parameter |1 -\- re^^e^'^^°"'\ ~ 1.9, since ^^ ~ 3.6 according to our 
FEM calculations. 

We now use the reflection coefficient to construct the full wave packet driving 
the atom, assuming that we initially launch a gaussian wavepacket 0j(z, t) = 
(27rA2)i/4e[i(fco2-^ot)]eh^^^^^|^±^2L] ^^^^ ^ ^ _^^ at t = into the waveguide. The 
backward propagating wavepacket due to the reflection can be written as 

MZ,t) = re^-^(27rA^)^/^e'^°"eW^'°^+-»^))et- "'"'""'^4^"°'"'\ (18) 

where Z = z — a. Hence the total fleld given by 0t(-2, t) = [4>i{z, t) + (f)r{Z, t)] at t = is 

Mz, 0) = (27rA2) V4e«-o*))e[-^^^4^] (19) 



+ re^'^(27rA2)i/4e^2'^'"'e-*'^o*e[- 

We simplify this expression using one assumption, namely that the amplitude envelope 
factor e' * ^ in the second term can be simplifled to be equal to the amplitude 

envelope factor e^ * i of the flrst term. The rationale for this approximation 
is that we consider single-photon wavepackets with a bandwidth of a few gigahertz, 
corresponding to an envelope length of a few tens of centimeters. The offset 2a in 
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the envelope factor of the second term, which is 2a=380 nm, is neghgible compared to 
the entire envelope length. In other words, the femtosecond travel time between the 
emitter and the wire termination is far below the nanosecond temporal length of the 
pulse envelope. As a consequence, we can take the envelope of the reflected pulse to be 
identical to that of the incident wavepacket. In this approximation, the total field that 
appears in the excitation term is 

0i(z,O)~0,(;2,O)[l + re^'(^+2'=°")]. (20) 

Exactly as in the case of a symmetric waveguide, the dynamics of the atomic excitation 
in Eq. IQ is set by the field amplitude 0i(fc,O), i.e., the Fourier transform of 0i(z,O), 
and the coupling strength g^ for the corresponding mode. Due to the reflection 
from the rod end, the field is enhanced by a factor of [1 + re-^'-^"'"^'^''"^], as evident 
in Eq. ( I20l) . Not only the field, also the coupling strength g^ is affected by the 
presence of the wire termination, since the decay rate into the plasmonic channel, 
i.e., Tpi = /^qF = 27r[l + re^^^~^'^''°"'^]'^gl/vg, again contains the reflection at the wire 
termination. The overall effect of the strong reflection at the wire end on the integral 
kernel in the excitation term is given by 

Cl^{0)gk = [1 + re^('+^''-%{k, 0)^ /3oTvg/{2TT[l + re^ie+2koa^2) (21) 

= 0,(fc,O)^/3ort;,/(27r), 

where 0j(A;,O) is the Fourier transform of initially incident wavepacket 0j(z,O). 
Equation fl?T]) shows that the factor [1 + re-^'-^"''^^'^"-'] cancels out in the excitation term. 
Therefore, the equation defined in Eq. (fTTj) has properly taken the reflection coefficient 
as well as the reflection phase into account, except that we assume the envelope to 
be constant over the spatial range from the emitter to the rod end, which is a valid 
approximation as discussed. 

4-3. Enhanced single-photon absorption by a terminated metallic nanowire 

We calculate the time dependence of the atomic excitation probability of an atom 
coupled optimally to the terminated metallic nanowire for a gaussian input photon 
wavepacket, as shown in Fig. [5] (c). For such a terminated metallic nanowire, we flnd 
that the maximum of the atomic excitation probability is 72.4%. The fact that the 
probability of atomic excitation exceeds 50% despite the fact that light is only injected 
into the system from one side, is a direct consequence of the high spontaneous emission 
beta factor of 91% for unidirectional emission into the open end of the waveguide. For 
an ideal case with /3-factor of 100%, we also flnd that the maximum atomic excitation 
is ~ 80% for a gaussian wavepacket. The fact that there is still a discrepancy between 
72.4% and 91%, (resp. 80% and 100%) shows that gaussian input wavepackets do not 
form the optimum temporal pules shapes. By pulse shaping the input wavepacket it 
may be possible to more perfectly approach the perfect inverting pulse. 

As an example of further optimization of pulse shape, we study the possibility 
of shaping the input wavepacket through waveguide dispersion. Considering a 
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highly dispersive waveguide with dispersion relation a{u) = ao + ai{ijJo){ijj — 
u) + a2{ujo)/2{u — uY + a2,{uJo)/6{uj — w)^, the gaussian wavepacket (f)i{—Zo,t) = 

(27rA^)^/'^ef~*'^°*]e[ * ^ initially launched at z = —Zq will be modified due to 

the waveguide dispersion. When the wavepacket reaches the emitter, the corresponding 
pulse shape is given by [17] . 

0,(0, t) = 4^(^)3/4L?li:^g[(2-3AB-6C^)^-/C(3AB+2C^-6)^] (22) 

r 
X Ai[{l -AB-C^ + I2C){\B\)-^/% 

where A = 4(t - aoZo)/r, B = 32Zoff , C = 8Zof|, r = 4:/{Avg), Ai{x) is the Airy 
function. By taking into account the dispersive features of the waveguide as well as 
the possibility of chirping the input pulse, we find that the chirping as well as a2 are 
simply detrimental to reaching maximum absorption, since they mainly give rise to the 
broadening of the pulse. Interestingly, an optimized a^luo) can improve the maximum 
of the atomic excitation by 4% compared with the nondispersive case, as shown in the 
contour plot of maximized atomic excitation as function of bandwidth (Q) and a^ in 
Fig. [5] (d). In closing, we have shown that considerable amount of atomic excitation can 
be gained simply by operating at the termination of the plasmonic waveguide. 

5. Conclusion 

In conclusion, we have theoretically analyzed the time-dependence of single-photon 
absorption by a single emitter. The model allows us to specify any single-photon 
wavepacket as initial condition, and to calculate the time dependence of the atomic 
excitation. We apply the theoretical model to quasi- ID waveguides coupled to a single 
emitter. For a gaussian single-photon wavepacket, we give the analytical solution to the 
atomic excitation, as well as the numerical results for transmitted and refiected light. 
We compare our time-dependent theory to the stationary theory, and we have reported 
that our time dependent theory contains the stationary refiection and transmission 
spectra of earlier work by Shen and Fan. To optimize the excitation probability of the 
emitter it is essential to choose incoming wavepackets of optimum duration. Within 
the class of gaussian single-photon wavepackets, excitation efficiencies up to 40% are 
possible for emitters coupled to infinite waveguides. We further studied the impact of 
the finite bandwidth of the incoming wavepacket on the atomic excitation, and find 
the maximum excitation probablity 40%. This high atomic excitation simply generated 
by the gaussian distributed wavepacket is close to the fundamental limit of 50% set by 
reciprocity. In order to obtain an even higher excitation probability, we have proposed to 
engineer the photonic environment of the emitter to suppress the spontaneous emission 
into all channels, except the one into which the incoming single photon is funneled. 
Practically, by terminating a plasmonic nanowire and positioning the emitter properly 
we have found that most of the light, up to 91%, can be directed into can be directed into 
a single channel. Reciprocity guarantees that if we use the high spontaneous emission 
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/3-factor channel for pumping, a high atomic excitation in excess of 50% can be achieved. 
Indeed, we find that a value of 72.4% can be achieved by using a very simple structure, 
i.e., a terminated metallic nanowire. This result is obtained with a gaussian distributed 
single-photon wavepacket, and can be improved by shaping the optical pulse. Using 
waveguide dispersion for pulse shaping, a further modest improvement can be achieved. 
As an outlook, we envisage that our time dependent-theory for the atomic excitation 
can be useful for the analysis of different experimentally relevant scenarios, and on- 
going experimental activities that focus on coupling a freely propagating photon to 
an atom fiU\ [TT| [12] . The relevance of engineering the photonic environment of the 
emitter is evident from our previous discussions for achieving high atomic excitation 
rate, and might provide new guidelines to do future experiments. Particularly promising 
for achieving high atomic excitation is to use 3D photonic structures, i.e., optical 
nanoantennas [HI |26l [27], where the temporal and spatial coherence of the incoming 
wavepacket can be fully addressed. Optimizing single-photon absorption using the 
methods discussed in this paper will have high impact throughout the photonics 
community, spanning from quantum optics, single molecule absorption microscopy, to 
photovoltaics. 
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Appendices 



Appendix A. Normalization in continuous-mode quantization scheme 



In the well-known discrete-mode quantization scheme j43|, the commutation rules for the 
photon operator obey [aA,ay] = 6xx'- The normalization condition for these discrete 
modes is J e{r)[Ex{r)] ■ [Ex'{r)]dV = Si3^,3>6p^p>6m,m'Sq,q'Ndis. Accordingly, the single- 
photon field in the discrete- mode quantization scheme is obtained as Ex,i,dis{f) = 
Ex{r)^/ ^^ , where u = a/c^(/3^ + q'^)- In this paper we deal with ID waveguides, 
and hence a quantization box with infinite extent parallel to the z axis but with a finite 
cross sectional area in the transverse plane. We extend Eq. ([6]) to the continuous mode 
quantization scheme. The mode spacing along the z direction, A/3, tends to zero as the 
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quantization length along z axis approaches infinity. Therefore, the sum over all the 
modes ^ can be substituted by 

A 

A {m,p,/3,<?} ' •' n ' •> ^^ > k 

where A/3 is the mode spacings along the z direction. In this substitution rule, 
^ represents the sum over all the modes with frequency u. The normalization 

K 

condition in the continuous-mode quantization scheme needs to be modified to read 
^ e{r)[Ex{r)]-[Ex'{r)]dV = S{l3 - l3')6p^p>6m,m'5g,g'Ncon- Accordingly, the single- 
photon field in the continuous- mode quantization scheme is obtained as Exicon{f) = 



E\{r)J ^^^ . As a side remark, the discrete Kronecker delta and the continuous Dirac 
delta- funtion are related by df^^pi — )■ A/3(5(/3 — /3'), which indicates that the normalization 
factor has incorporated a factor of A/3 that yields the translation of the single-photon 
field from the discrete mode quantization to the continuous mode quantization as 

Ex,i,con{r) = Ex,iMr)/V^- 

Appendix B. Laplace transform method for gaussian input single-photon 
wavepacket 

An alternative approach to obtain the atomic excitation probability as a function of 
time that is distinct from our direct time integration of Eq. (j8]) to obtain Eq. ( TT2|) is 
to use a Laplace transform method reported by Wubs et al. [44J. In this appendix we 
summarize the Laplace method. By Laplace transforming the Heisenberg's equation of 
motions of the atomic operator h{t) and field operator a(t), one can include the initial 
atomic excitations h{t = 0) and photonic excitations a{t = 0) in the dynamics. More 
details can be found in [44]. After making a pole approximation, the inverse Laplace 
transform gives the dynamical solution to the atomic operator, 

b{t) = 5(o)e-^(^^+^)*-r*/2 (B.l) 

1 r-oo -iuj't _ -i{nA+e)t-Tt/2 



1 f~°^ 



where F is the total spontaneous-emission decay rate, e is the Lamb shift of atomic 
transition Qa of the emitter due to the coupling to the waveguide modes, g{uj') is 
the coupling strength, 6(0) is the initial atomic excitation, and a(t^/^+)(0) (a((^/^+)(0)) is 
the initial right (left) propagating optical excitation. The atomic population operator 
N^it) is given as (A^a(^)) = {i\b'^ (t)b(t)\i) , with the initial condition given by \i) = 
\g) (8) J dijj'S{uj')d(^i^r^+){0)\0)p, and S{uj') defines the incident single-photon wavepacket. 
Using the same gaussian distributed single-photon wavepacket as initial condition that 
is also used in Sec. 3, i.e., S{k) = Cl{t = 0), one finds \i) = \g) ® j dkS{k)a\ _^\Q) , 
which can be reformulated as an integral over frequency \i) = \g) ® J duj'xi^')(\jt _^|0), 
with xi^) = S{^)/ \/'^g{'-^)) (4u = a{./A/fg(w). The group velocity and phase velocity 
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are defined as Vg = du/dk, Vp = u/k. Within the linear dispersion approximation, Vg is 
equal to Vp. The expectation value of the excited atomic state reads as follows, 



{NA{t)) 



dug, 



^-iujt _ g- 



-i(uiA+e)t~ir)t/2 



-x(w) 



Lu-nA-e + i{T)/2 
Using a = u — Uo and omitting the Lamb shift term e, one obtains 

I d.rr 

47r 



{NA{t)) 



TX 



fi2 



'^^^^/"" a + .r/2 




(B.2) 



(B.3) 



By rewriting ]p,^ as ^ Jp,^ = —i J^ d^e''^^'^^^^'^^^ and performing the a-integral, one 
can further simplify (A^^(t)) as follows. 



{NA{t)) = K{t)\' = 
which is consistent with Eq. ( IT2l) . 



4A/^e#-[erf(v/^t + A) _ erf( A)]e-i* 
2 y a A/a i/a 



, (B.4) 
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